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ABSTRACT 

We test semi- analytic models for galaxy formation with accurate kinematic data of damped Lya pro- 
togalaxies presented in a companion paper (Wolfe & Prochaska 2000). The models envisage centrifugally 
supported exponential disks at the centers of dark matter halos which are filled with ionized gas undergo- 
ing radial infall to the disks. The halo masses are drawn from cross-section weighted mass distributions 
predicted by CDM cosmogonies, or by the null hypothesis that the dark matter mass distribution has not 
evolved since z ~ 3 (i.e., the TF models). In our models C IV absorption lines detected in damped Lya 
protogalaxies arise in infalling ionized clouds while the low ion absorption lines arise from neutral gas in 
the disks. Using Monte Carlo methods we find (a) The CDM models are incompatible with the low ion 
statistics at more than 99% confidence whereas some TF models cannot be ruled out at more than 88% 
confidence, (b) Both CDM and TF models are in general agreement with the observed distribution of C 
IV velocity widths, (c) The CDM models generate differences between the mean velocities of C IV and 
low ion profiles that are compatible with the data, while the TF model produces differences in the means 
that are too large, (d) Both CDM and TF models produce ratios of C IV to low ion velocity widths that 
are too large, (e) Neither CDM nor TF models generate C IV versus low ion cross-correlation functions 
compatible with the data. 

While it is possible to select model parameters resulting in agreement between the models and the 
data, the fundamental problem is that the disk-halo configuration assumed in both cosmogonies does not 
produce significant overlap in velocity space between C IV and low ion velocity profiles. We conjecture 
that including angular momentum of the infalling clouds will increase the overlap between C IV and low 
ion profiles. 

Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 



1. INTRODUCTION 

This is the second of two papers which discuss ion- 
ized gas in damped Lya systems. In Paper I (Wolfe & 
Prochaska 2000) we presented velocity profiles drawn from 
a sample of 35 damped Lya systems for the high ions C 
IV and Si IV and the intermediate ion Al III. Comparison 
among these profiles and with profiles previously obtained 
for low ions such as Fe II showed the damped Lya sys- 
tems to consist of two distinct kinematic subsystems: a 
low ion subsystem composed of low and intermediate ions 
and a high ion subsystem containing only high ions. The 
evidence distinguishing between the kinematic subsystems 
is robust and stems from a battery of tests comparing dis- 
tributions of various test statistics. It also stems from 
differences between the C IV versus low ion or C IV ver- 
sus Al III cross-correlation functions on the one hand, and 
the C IV versus Si IV or Al III versus low ion cross correla- 
tion functions on the other. Whereas the latter have high 
amplitude and small velocity width, the former have lower 

x Visiting Astronomer, W.M. Keck Telescope. The Keck Obser- 
vatory is a joint facility of the University of California and the Cali- 
fornia Institute of Technology. 



amplitude and wider velocity width. This is because ve- 
locity profiles of ions arising in the same kinematic subsys- 
tem comprise narrow velocity components that line up in 
velocity space, whereas velocity components arising from 
ions in different kinematic subsystems are misaligned in 
velocity space. However, the existence of a statistically 
significant C IV versus low ion cross-correlation function 
suggests the two subsystems are interrelated. This is also 
indicated by a systematic effect in which the C IV profile 
velocity widths Awciv are greater than or equal to the low 
ion profile velocity widths Aui ow in 29 out of 32 systems. 
In Paper I we claimed these phenomena indicate that 
the two subsystems are located in the same gravitational 
potential well. In this paper we shall expand on this idea 
with specific models. The models assume that the low 
ion subsystems are centrifugally supported disks of neu- 
tral gas located at the centers of dark matter halos (see 
Mo et al. 1998; hereafter referred to as MMW), whereas 
the high ion subsystems comprise photoionized clouds un- 
dergoing infall from a gaseous halo to the disk. That is, 
we assume the dark matter halos contain hot virialized 
gas in pressure equilibrium with the photoionized clouds. 
The hot gas undergoes a subsonic cooling flow toward the 
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disk while the denser clouds infall at velocities approach- 
ing free fall (Mo & Miralda-Escude 1996; hereafter referred 
to as MM). The models are set in a cosmological context 
by computing the mass distribution and other properties 
of the dark matter halos using Press-Schecter theory and 
the CDM cosmogonies adopted by semi-analytic models 
for galaxy formation. We also consider the null hypothesis 
that galaxies at z w 3 have the same dynamical properties 
of nearby galaxies (this model, hereafter referred to as the 
TF model, is defined in § 2.4.2). 

We test the models using a Monte Carlo technique for 
computing absorption spectra arising when sightlincs that 
randomly penetrate intervening disks also intercept ion- 
ized clouds in the halo. § 2 presents models for the neutral 
gas. We discuss properties of the disk, the cosmological 
framework, and Monte Carlo techniques. In § 3 we dis- 
cuss models for the ionized gas. We derive expressions for 
the structure and kinematics of the two-phase halo gas. 
Expressions for C IV column densities of the clouds are 
derived. In § 4 we give results of the Monte Carlo simu- 
lations for the low ion gas for both CDM and TF model. 
Results of Monte Carlo simulations for the ionized gas are 
given in § 5. Here we also consider tests of correlations 
between the kinematics of the high ions and low ions. In 
§ 6 we investigate how the results of § 5 are affected by 
changes in some key parameters such as central column 
density and low end cutoff to the input dark-matter halo 
mass distribution. Concluding remarks are given in § 7. 

2. MODELS OF THE NEUTRAL GAS 
2.1. Cosmological Framework 

To place the model in a cosmological context we as- 
sume bound dark-matter halos evolve from linear density 
contrasts, 6(x,t)=6p/p, according to gravitational insta- 
bility theory for Friedmann cosmologies (Peebles 1980). 
We consider adiabatic CDM models (ACDM) in which 6k , 
the Fourier transforms of 6(x, t), are Gaussian distributed 
with variance given by P{k), the power spectrum at the 
epoch of radiation and matter equality. We also consider 
isocurvature CDM models (ICDM) in which the 6k are not 
Gaussian distributed (Peebles 1999b). The 6k, or more 
specifically the rms density contrasts in spheres with mass 
scale M, i.e., Am, grow with time until they go non-linear 
and collapse. According to the spherical collapse model, 
this occurs when the Am predicted by linear theory equals 
6 C = 1.68. To compute n(M,z)dM, the density of bound 
halos in the mass interval (M, M+dM), we follow pre- 
vious authors (e.g. MMW) who used the Press-Schecter 
formalism in the case of ACDM. In Appendix A we derive 
an expression for n(M, z)dM in the case of ICDM. 

We shall also test the null hypothesis that little or no 
evolution of galaxies has occurred since the epochs of damped 
Lya absorption; that is, an hypothesis assuming current 
disks, with higher ratios of gas to stars, were in place at 
z > 3. This scenario resembles the semi- analytic models 
in that we assume centrifugally supported disks reside at 
the centers of dark matter halos filled with hot gas at the 
virial temperatures. However, we do not assume a CDM 
power spectrum nor the Press-Schecter formalism to com- 
pute the mass distribution of halos. Rather, we assume (a) 
the luminosity function of galaxies in the redshift range of 
the damped Lya sample in paper I, i.e., z — [2,5], is given 



by the Schecter function determined from nearby galaxies 
(e.g. Loveday et al. 1999), (b) galaxies at these redshifts 
obey the same Tully-Fisher relationship between luminos- 
ity and rotation speed as at z — 0, and (c) a correlation 
between disk radial scale-length and disk rotation speed 
exists. This model, hereafter referred to as TF, is an ex- 
tension of the rapidly rotating disk model suggested by 
Prochaska & Wolfe (1997 and 1998; hereafter PW1 and 
PW2). We describe this model in more detail in the fol- 
lowing sections 2 

Throughout the paper we shall intercompare results for 
the TF model and four CDM cosmogonies. The cosmo- 
logical settings of the CDM models are specified by (i) 
the current total matter density, Am, (ii) the cosmologi- 
cal constant, fl\, (iii) the Hubble constant, h, where h = 
Ho/ 100 km s _1 Mpc -1 , (iv) as, the rms linear density 
contrast at z = in spheres of radius 8ft- -1 Mpc, and (v) 
n, the power-law index for the power spectrum, in cases 
where P(k) ex k n . The values of the parameters are given 
in Table 1. The SCDM, ACDM, and OCDM are normal 
ACDM models. In all three cases, P(k) is given by the 
Bardeen et al. (1986) expression (which is not a power 
law). In the ICDM model we follow Peebles by assuming 
P(k) oc k n and n = —1.8. The cosmological parameters 
specifying the TF model are also given in Table 1. 

2.2. Disk Models 

Most semi-analytic models assume the neutral gas caus- 
ing damped Lya absorption is confined to centrifugally 
supported disks at the centers of dark-matter halos (e.g., 
Kauffmann 1996; MMW). The spherical collapse model is 
used to define the limiting virial radius as r2oa, the ra- 
dius within which the mean density of dark matter equals 
200p cr it(z) where p cr it{z) is the critical density of the uni- 
verse at redshift z. Thus, V200, the circular velocity at 
?"20Cb is related to r2oo and M, the halo mass within r2oo, 

by 



V 2m = [10GH(z)M] 1 / 3 ; r 200 = 



V: 



200 



10H(z) 



(1) 



where H(z) is the Hubble parameter, (a~ 1 da/dt), at z. 
MMW also define rrid and jd as the fractions of halo mass, 
M, and angular momentum, J, in disk baryons. Assum- 
ing the halos to be singular isothermal spheres embedding 
disks having exponential surface-density profiles with ra- 
dial scale- lengths, Rd lso , they find that 




(2) 



and 



2 We could have adopted luminosity functions that arc measured 
in this redshift interval. The most accurate determinations are for 
the Lyman- break galaxies (Stcidcl ct al. 1999). However, we re- 
jected this procedure because the rotation speeds of these objects 
have not been measured, nor is it known whether or not they con- 
tain rotating disks. The strong clustering exhibited by the Lyman- 
brcak galaxies suggests otherwise. Nevertheless, we briefly discuss 
this possibility in § 6.2 
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N^ so = 4.8x10 




I 250kms" 1 



(3) 



where the spin parameter of the halo, A = J\E\ 1 / 2 G^ 1 M^ b / 2 , 
(E is the total energy of the halo) and Nq is the central H 
I column density perpendicular to the disk. The distribu- 
tion of A is given by 



p{X)d\ = 



1 



2it<t\ 



exp 



ln 2 (A/ < A >) 



2a\ 



f , (4) 



where the mean and dispersion are determined from nu- 
merical simulations to be < A > = 0.05 and o\ =0.5 (see 
Barnes & Efstathiou 1987). 

For the TF models we infer the parameters of the halo 
from observed properties of the disk. Thus, we use model 
rotation curves to infer V200 from V max , the observed max- 
imum rotation speed (see below). We then use eq. (1) 
to obtain the mass and virial radius of the halo. In this 
case we obtain the radial scalelength and central column 
density of the disk by adopting the following correlations 
inferred by MMW from the spiral galaxy data of Courteau 
(1996; 1997): 



R 



Cor 



-1.25 + 7.4(y maa; /250kms _:L ) hT 1 kpc 



(5) 



"'V -1 h cm" 2 , (6) 
where [i is the mean molecular weight of the gas. 



N Cor = 3 6xl0 22+0 - 45(Vmax / 250kn 



2.3. Rotation Curves 

We next turn to the mass distribution of the halo. This 
is crucial for determining both the rotation curve of the 
disk and the dynamics of gaseous infall discussed in § 3. 
For our model we adopt the analytical fit to the halo mass 
distribution found in the N-body simulations of Navarro et 
al. (1997; hereafter NFW). In this case the halo rotation 
curve, defined by V ro t{r)=^/GM{^r)/r 1 has the following 
form: 



V rot (r) = V, 



2001/ — - 



c ln(l + x) - x/(l + x) 



r/r s 



where the concentration parameter, 



x m(l + c) - c/(l + c) ' 

(7) 
7*200 /r s , and at 
r = r s , the halo mass density p ex r~' 2 . NFW developed 
a self-consistent theory in which c depends on V200 and 
redshift, z, in the sense that at a given z, c declines with 
increasing V200, and at a given V200, c decreases with z. 
We use the algorithms described in Navarro et al. (1997) 
to compute c=c(z, V2oo)- NFW halo rotation curves corre- 
sponding to z = 2.5, A = 0.05, m = 0.05, and the ACDM 
cosmology appear as solid curves in Figure 1 for V200 = 
50—250 km s _1 . As expected V ro t — V ma x at r w 2r s . 
Curves with lower V200 appear to rise more rapidly in the 
interval r = [0,r2oo] because the ratio r s /r2oo decreases 
with decreasing T^oo- 



However, the expression for the halo V ro t(r) in eq. (7) 
is incomplete as it ignores the presence of the disk. Self- 
gravity of the disk affects the mass distribution of the halo 
through adiabatic contraction (Blumcnthal et al. 1984). 
We used the MMW formalism to compute rotation curves 
due to contracted halos and found that V rot (r) differed 
from the expression in eq. (7) by less than w 15%. Given 
the uncertainties in the models we shall ignore these cor- 
rections and use the halo mass distribution implied by eq. 
(7) when computing dynamics of infall . 

On the other hand the rotation curve of the disk can 
differ significantly from eq. (7) when the disk contribu- 
tion to the net potential gradient is added to that of the 
contracted halo. MMW compute the scale length and cen- 
tral column density of centrifugally supported exponential 
disks in adiabatically contracted NFW halos formed by 
spherical collapse. In comparison with halos modeled as 
singular isothermal spheres they find 



pMMW 



-V2, 



-RTf c ' irz fR ,N MMW =Nrfcfu Z ■ (8) 
Explicit expressions for the functions f c and /# are given 
by MMW. They then use the new disk parameters to cal- 
culate the modified V ro t(r) for disks embedded in NFW 
halos. Examples of modified V ro t(r) corresponding to A = 
0.05, md = 0.05, and z — 2.5 are shown as dotted curves in 
Figure 1. Clearly the rotation speeds sampled by sightlincs 
traversing these model disks lie between l^oo and iyV^ooi 
the maximum of the modified disk rotation curve, where 
fv, a function tabulated by MMW, exceeds unity. Given 
the wide range of possible rotation curves appropriate for 
model damped Lya systems (see PW2), we assume that 
the disk rotation curves have either of the following con- 
stant speeds: 



V ro t(r) = V200X 



1 

fv 



(9) 



For the TF models V200 is not given apriori . Rather we 
assume V ro t (r) = V max and that V max is selected from an 
empirical distribution derived from the Tully-Fisher re- 
lation (see § 2.4.2). The radial scale length and central 
column densities then follow from eqs. (5) and (6). In this 
case V200 = Vmax or V200 =V max /f v . 

2.4. Monte Carlo Models 

In previous work (PW1, PW2) we tested models by com- 
paring predicted and empirical distributions of the four 
test statistics for the low ions (Figure 6 in Paper I). The 
model distributions were computed by a Monte Carlo tech- 
nique in which low ion absorption profiles were produced 
by sightlines traversing 10000 randomly oriented disks, 
and test statistics were determined for each profile. In 
PW1 the disk models were based on the simplifying as- 
sumption of identical disks with flat rotation curves char- 
acterized by a single rotation speed, while in PW2 more 
realistic forms of V ro t{r) were used for the identical disks. 
Here we extend this approach to account for the distribu- 
tions of halo masses and spin parameters. 

2.4.1. CDM 

MMW give an expression for the cumulative probability 
that sightlincs to the background QSOs intercept disks in 
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Radius (kpc) 

Fig. 1. — NFW rotation curves corresponding to z = 2.5, A =0.05, 
m = 0.05 and the ACDM cosmology. Solid curves are halo rota- 
tion curves from eq. (9) corresponding to V200 = 50, 100, 150, 
200, and 150 km s _1 . Dotted curves show rotation curves re- 
sulting from modifications due to adiabatic contraction and disk 
self-gravity. Curves plotted out to r^oo 



halos with masses exceeding that corresponding to V max '- 
the result is averaged over p(A). For the Monte Carlo 
model we require the differential expression; that is, the 
probability for intercepting disks in the spin parameter 
interval (A, A + dX), circular velocity interval (V200, V200 + 
rfV^oo)) and redshift interval (z, z + dz). The result is given 

by 



H I column densities less than Ni are thrown out. We re- 
strict the boundaries of the surface to V200 < 300 km s _1 
to insure that gas in virializcd halos has ample time to cool 
and collapse to the disk (e.g. Rees & Ostriker 1977) by 
a redshift, z = 2.6, the median redshift of the kinematic 
sample. We also assume dPcDM cuts off below V200 = 30 
km s _1 since gas photoionized to temperatures of ~ 10 4 
K by the UV background radiation escapes from the dark 
matter halos with V200 < 30 km s _1 (Thoul & Weinberg 
1995; Navarro & Steinmetz 1997; Kepner et al. 1997: we 
investigate the consequences of modifying this restriction 
in § 6). Figure 2 shows the resulting distributions of V200 
for all the CDM models in Table 1. All the CDM curves 
exhibit maxima near the V200 cutoff predicted by hierarchi- 
cal cosmologies, and as predicted the curves decline with 
increasing V200 • As expected the largest fraction of mas- 
sive halos is indicated for the ICDM models. Thereafter, 
the fraction decreases progressively from OCDM, ACDM, 
to SCDM adiabatic models. 



Table 1 

Model Parameters 

MODEL SCDM ACDM OCDM ICDM TF 



S2m 


1.0 


0.3 


0.2 


0.2 


0.3 


n A 


0.0 


0.7 


0.0 


0.8 


0.7 


h 


0.5 


0.7 


0.7 


0.7 


0.7 


C"8 


0.6 


1.0 


1.03 


0.9 




n 


a 


a 


a 


-1.8 





"CDM power spectrum given by Bardeen et al. (1986) expression 



2.4.2. TF 



dP C DM{z,\,V 2 oo)= dz|(l + z) 3 (-c— ) 



xdXp(X) [RT(V 2 oo, A, z)N™° (V 200 , A, z) 

xdV wo n{V 2Q o,z)F(N^ ) , (10) 



where 



F(N l so ) = 



and 



{Nr)2 



\ + HN^/Nt) (l + ln(JV2"°/JV,) 

(11) 



n(y 20 o,z) = n(M,z) 



dM 



dV 2( 



(12) 



Here Ni = 2x10 cm is the threshold column density 
for damped Lya surveys (e.g. Wolfe et al. 1995). 

For given redshift intervals and cosmologies, the func- 
tion dPcDM describes a surface above the (A, V200) plane. 
To form synthetic samples of 10000 low ion profiles, we 
randomly draw (A, V200) pairs according to the height of 
the surface above the plane. To assure compliance with 
the damped Lya surveys, sightlines resulting in observed 
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Fig. 2. — Input circular velocity distributions for the CDM and TF 
models used in Monte Carlo simulations. Computed as described in 
the text for the velocity interval (30, 300) km s — 1 . The abscissa 
correspond to V200 in CDM and Vmax in TF models. 
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The curve labeled TF in Figure 2 represents the null 
hypothesis discussed above. In this case the x-axis cor- 
responds to Vmax rather than V-zoo- This is because V200 
is a theoretical construct, whereas the null hypothesis is 
based only on the observed properties of current galaxies. 
The crucial relation ship here is the Tully-Fisher equation 
which connects L and V max ; i.e., 



L(V n . 



L*{v max /v*y 



(13) 

where 14 is V max for a galaxy with luminosity L — L*. 
With the last equation we can obtain expressions for the 
interception probability 



dt. 



.^2 (1+Z)3 (- C i) 



dP T F(z,V max ) 

x l^d ("mon 2)-/Vq (V max , Z) 



XdVrr 



n(V max ,z)F(N^° r ) , (14) 

where dV max n[L(V max )] is the present density of galaxies 
with maximum rotation velocities in the interval (V max , 
V max +dV max ), and F(Nfr' or ) is given by eq. (11) with 
N§ or substituted for N^ so . To determine n[L{V max )} we 
assume a Schectcr luminosity function; i.e., 



dV max n(V max , z) = $»(£/£») a cxp(-L/L*)d(L/L*) . 

(15) 
The TF curve depends on the parameters a, 7, and 14 
for which we assumed values of 1.0, 3.0, and 250 km s _1 
respectively. These are representative for the values of 
parameters adopted by Gonzalez et al. (2000) who com- 
puted n(V max ,0). They analyzed data sets from exten- 
sive surveys carried out at B magnitudes. However, using 
Ccpheid calibrated galaxies Sakai et al. (2000) derive 
Tully Fisher relations indicating lower values of 14 (« 180 
km s _1 ) for the B magnitude Tully Fisher relation. Low 
values of 14, as well as higher values of 7, are indicated by 
their / band Tully Fisher relation, and from similar rela- 
tions found by Giovanelli et al. (1997). In § 6 we discuss 
the sensitivity of our results to these parameters . The 
point we wish to emphasize is that for acceptable ranges 
of these parameters the TF curve differs from the CDM 
curves in that it (a) peaks at V max « 14/2 s=s 130 km s _1 
, which is large compared to the 30 km s _1 peak of the 
CDM curves, (b) has little power near the latter peak, and 
(c) falls off exponentially when V200 > 14/2. Of course the 
comparison with CDM is inexact since Pcdm depends on 
V200 while Ptf depends on V max . Still the differences be- 
tween V200 and Vmax are not large enough to invalidate 
these conclusions. 

3. IONIZED GAS 

Assume that disks arise from the infall of ionized gas 
predicted to fill dark-matter halos. According to the MM 
model the halo gas is accreted during merger events with 
other halos and is shock-heated to the virial temperature 
of the halo, kT V i r = (l/2)(/zmjyV2 2 oo); we assume [i = 0.4. 
The duration of the accretion phase is presumably short 
compared to tM, the time- interval between events. For the 
mass interval corresponding to V200 = [30, 300] km s _1 , 
T V i r ranges between 3xl0 4 and 3xl0 6 K. In the case of 
massive halos the cooling time, t coo i, exceeds the age of 



the gas, tM, at r — r2oo- Because the density of the gas 
increases with decreasing radius, t coo i decreases with de- 
creasing radius until t coo i = tM at the cooling radius, r coo /. 
At r < r coo i the gas moves radially inward in a quasistatic 
cooling flow (Fabian 1994). Cool clouds form in pressure 
equilibrium with the hot gas since the hot gas is thermally 
unstable. Due to the loss of buoyancy, the denser clouds 
fall inward at speeds determined by the imbalance between 
gravity and ram pressure. Because of rapid cooling, most 
of the gas in lower mass halos cools before moving inward, 
but within a limited range of halo masses there is always 
hot low-density gas left over which also moves inward in 
a cooling flow and exerts pressure on cooler clouds which 
have formed. The MM hypothesis is that the cool clouds 
are photoionized by background UV radiation and that 
they are the sites of C IV absorption lines in QSO absorp- 
tion systems. Here we extend this hypothesis to model 
the ionized gas causing C IV absorption in damped Lya 
systems. 
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Fig. 3. — Solid curves are pressure profiles for hot gas in halos 
with V200 = 100 and 250 km s _1 in (a) and (b) respectively. Pres- 
sure profiles are given by product of p(r) h and 77 (r), the density 
and temperat ure of the hot gas. Profiles computed according to 
prescription in Appendix B. Dot-dashed curves are smooth density 
of cool gas times temperature of clouds; i.e. p c T c i d . When latter 
exceeds pressure of hot gas, prcsurc equilibrium between cool clouds 
and hot gas breaks down (see § 4.1) 



3.1. Two Phase Structure of the Halo Gas 

Following MM we assume the hot gas is in hydrostatic 
equilibrium with the dark matter potential and that it ex- 
hibits adiabatic temperature and density profiles out to 
r m in = min(r coo/ ,r 20 o) • At r > r mrn the hot gas either 
follows isothermal profiles when r2oo > f coo i or does not 



exist when r2oo < 



In Appendix B we derive expres- 



sions for r coo i, and for ph(r) and T n (r), the density and 
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temperature profiles for hot gas in halos corresponding to 
the unmodified NFW rotation curves in eq. (7) . Examples 
of pressure profiles for high and low mass halos are shown 
in Figure 3 . The Figure shows the gas pressure in high 
mass halos to be 10 to 100 times higher than in low mass 
halos. The high pressures are mainly due to the higher 
virial temperatures of massive halos. 

Gas in the cool phase comprises identical uniform spher- 
ical clouds with mass, M c u, radius, R c id, temperature, 
T c id, and internal density, pi nt - We adopt the MM model 
by assuming clouds form at r — r m i n with mass M c id(r m i n ) 
= 4xl0 5 Mq, and temperature, T c id =2xl0 4 K, since they 
are assumed to be photoionized. To compute R c id, MM as- 
sume pi n t is set by pressure equilibrium with the surround- 
ing hot gas; i.e., pint = PhThl^cM- But this assumption is 
incorrect for a wide range of halo masses and radii. Specif- 
ically, pressure equilibrium breaks down when the pressure 
is sufficiently low for R c id to equal the mean distance sep- 
arating clouds. This will occur when p c (r), the average 
density of cool gas, exceeds Pi n t{r)', i.e., when p c T c id > 
PhXh- According to MM, p c {r) is given by 



< \ ^1 



(16) 



where the constant C\ is evaluated in Appendix C. Fig- 
ure 3 demonstrates how r crossl the radius at which p c T c id 
= PhTh, decreases with increasing V200 In Figure 4 we plot 
r cross, r coo i, and r 2 r>o versus V" 2 oo for a ACDM cosmology, 
NFW halos, and z = 2.5 (we also plot the cooling radius 
of a singular isothermal sphere, R CO oi, for comparison with 
NFW halos). The point of the Figure is to show that (1) 
pressure equilibrium breaks down throughout halos with 
V200 < 100 km s _1 , (2) pressure equilibrium breaks down 
only near the centers of halos with V200 > 100 km s _1 , 

(3) r m i n = r2oo for halos with V200 < 250 km s _1 , and 

(4) r m i n — r coo i for halos with V200 > 250 km s _1 . These 
trends are qualitatively similar for all the background cos- 
mologies in Table 1. 

In order to compute the structure of the cool clouds we 



assume pi nt — p c at r < r c 



As a result 



Pint(r) = 



Ph(r)T h (r) 



:r > r c 



(17) 



Pc\' ) ' "^ ^cross 

The cross-section of the infalling clouds is then given by 



A(r) 



3M cW (r) 



2/3 



T c 



4tt 



y T h (r)p h (r) 
( I W3 



-p c {r)> 



)2/3 ■ r > r c 
, T -C T cross , 



which will be smaller than the value computed from pres- 



sure equilibrium at r < 



At r > r r 



M, 



id is 



assumed to be a function of r as the clouds may sweep up 
hot gas as they fall inward (see Benjamin & Danly 1997 
for a discussion of this problem) . At r < r cross the den- 
sity of hot gas is negligible, and so we assume M c id(r) = 
M c id(r cross ). Therefore we compute M cld {r) by solving 



dM c u = ( -A(r)p h (r) ;r > r c 
dr \ ;r < r cross , 



(19) 



subject to the boundary value of M c id(r cross ). Explicit 
expressions for A(r) and M c id(r) are given in Appendix 
C. 
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Fig. 4. — Critical radii versus V200 for ACDM cosmology, NFW 
halos, and z = 2.5. R coo i, the cooling radius of an infinite isothermal 
sphere is computed from eq. (B7), and r coo ; the cooling radius of an 
NFW halo is computed from eq. (B6). The "virial radius" r^oo is 
computed from eq. (3). The crossing radius r croBa , computed from 
eq. (C8), is the radius within which pressure equilibrium breaks 
down. Note r min = min(r 2 oo,?*cooi) for values of V200 where r cross 



< r 2 oo- At V200 < 100 km s" 



r-200 (see § 4.2). 



3.2. Cloud Kinematics 

The hot gas surrounding the clouds will exert a drag 
force opposing their radial infall. Assuming the drag is 
caused by momentum imparted to the clouds by hot gas 
swept up during infall, we find the radial equation of mo- 
tion to be 



d 

dt 



M dd V r 



= M cld g{r) 



(20) 



where g(r) is the gravitational acceleration and the radial 
velocity, V r = —dr/dt, is positive for infalling bodies. It 
follows that 



dV r 2 fdlnMj, 



dr 



{- 



1,1 



dr 



V 2 



-Mr) 



(21) 



(which corresponds to the case Cd = 2 of Benjamin & 
Danly [1997]). The solution to equation (21) is 



K 2 (r) = 



M 2 cld (r) 



drg{r)M 2 cld (r) + V?{r start ) 



(22) 

where we have computed the radial velocity of a cloud 
that starts to infall from an initial radius, r=r star t, with 
an initial velocity, V r (r sta rt)- Note, this expression ignores 
fragmentation of the clouds due to Kelvin-Helmholtz in- 
stabilities which may be important (Lin & Murray 2000). 
To compute V r (r) we assume the clouds infall from rest 
at r~r rn i n . As a result we solve eq. (22) with r sta rt — r m i n 



WOLFE & PROCHASKA 



0. The resultant V r (r) are valid in the 
Train]- The solution for V r (r cross ) acts as 



and V r (r start 

interval r — [r cross?' mm 

boundary condition in the first of two scenarios we consider 

for obtaining V r (r) at r < r cross . 

In the first scenario, clouds at r<r croS s follow pressure- 
free ballistic trajectories along which M c id(r) =const. In 
halos with r cross less than r m in we solve eq. (22) by assum- 
ing r start = r C ross and let V r {r start) equal the V r (r cross ) 
obtained from the solution at r > r cross . For halos with 
fcross greater than r m in we solve eq. (22) by assuming 
rstart = r m in and V r {r st art )=0; i.e., the clouds undergo 
ballistic infall from rest at r—r m in- 

In the second scenario we assume the cloud kinematics 
at r < r C ross are dominated by random motions. These 
may be generated by feedback from supernova remnants 
arising from star formation stimulated in cloud-cloud col- 
lisions. Such scenarios have been suggested to solve the 
"cooling catastrophe" characterizing the hierarchical build 
up of galaxy-scale structure in most CDM models (White 
& Frenk 1991). We obtain the velocity dispersion of the 
clouds by solving the Jeans equations for o~ r , the radial ve- 
locity dispersion for a system of clouds with (i) an isotropic 
velocity distribution, (ii) an average density distribution 
given by p c (r), and (hi) the gravitational field of an NFW 
halo (see eq. D4 in Appendix D). We then randomly draw 
the velocities of individual clouds from a Gaussian velocity 
distribution with dispersion given by ay. 

We emphasize that our model is most uncertain for low 
mass halos. This is because the underlying assumption of 
pressure equilibrium, which allowed MM to compute the 
properties of the clouds, breaks down throughout the in- 
falling gas for halos with V200 < 100 km s _1 . In these halos 
the properties of individual clouds are difficult to compute 
because without a confining medium the clouds become 
indistinguishable at r < r cr0 ss ■ Our approach to this prob- 
lem is to assume the line of sight traverses a medium con- 
taining "cloud-like" structures with fixed masses and that 
these give rise to C IV absorption lines. This assump- 
tion needs to tested with high-resolution hydrodynamical 
simulations of gas at r < r cr0 ss which is subject to input 
of mechanical energy. We have also assumed eq. (17) is 
valid to obtain the average density of the cool gas in ev- 
ery model. While this expression, which is based on mass 
conservation for infalling gas, is physically justified in the 
case of systematic infall, it is arbitrary when the gas kine- 
matics are dominated by random motions. Nevertheless 
we believe the results should provide insights into the ob- 
servational consequences of random motions (see § 6). 

3.3. Absorption Properties of the Cool Gas 

In order to compute C IV absorption profiles produced 
by the infalling clouds we select their locations along the 
line of sight from the cumulative interception probability 
function 



where the path integral propagates along a sightline with 
impact parameter b (where b is the distance in the plane 
of the sky separating the QSO sightline from the center 
of the galaxy), s max = v / r^ m - & 2 , and s mm = -s max . 
We compute N^f v (r), the C IV column density of a given 
cloud, from the following expression: 

N c c l f v (r) = X CIV (r)N c J d , (24) 

where the total H column density, i.e. H° + H + , of the 
cloud is given by 



N cld 



(r) 



3V2M c id(r)S 
AptmuAir) 



(25) 



and Xciv{ r )i the ratio of C +3 to total hydrogen volume 
densities, is assumed to depend only on r, the distance 
of the cloud from the center of the galaxy. To account 
for variations in N^f v caused by the variations of sight- 
line locations across the projected face of the cloud we 
have introduced the uniform deviate S which selects ran- 
dom numbers from the interval S =(0,v2). We evaluate 
Xciv( r ) by assuming 

X C iv(r) = X C iv(Rd)(Rd/rf . (26) 

To determine Xciv(Rd) and j3 we compute the average C 
IV column density, Nciv{b), where 



N C iv{b) 



We find that 



dsX CIV {r{s)) 



Pc(r(s)) 



pm H 

r = Vs 2 + b 2 . (27) 



where 



K = 



N C iv(b) = K\ 



2Ci X(R d ) 



R d 



1+/3 



dx 



(28) 



(29) 



pmH Rd Jo {1 + xZy+fJ/z ' 

and where x ma x = [(r m in/b) 2 — l] 1 / 2 . We then require 
the total C IV column density in n c id clouds to agree with 
Nciv{b) within some accuracy, e; i.e., . 



E 



n cld ATdd 

iV c/y 



in) 



Nciv(b) 



- 1 



< e 



(30) 



where e = 0.2 We adjust the input parameters K and /3 
by comparing empirical and model generated frequency 
distributions of C IV column densities (sec Figure 6) , and 
compute X(Rd) from eq. (29). 



F(< s) = 



r 


in dy 


Pc{r(y)) 
M cld {r{y)) 


A(r(y)) 


JSm, 


Z dy 


Pc(r(y)) 
M c u(r(y)) 


A(r(y)) 



= Vy T Tv 2 



(23) 



4. RESULTS OF MONTECARLO SIMULATIONS: LOW-ION 

GAS 

4.1. CDM 

Assume the low ion gas to be in disks with rotation 
curves normalized to V200 at r — ^200, and V200 to be ran- 
domly drawn from the distributions in Figure 2. In that 
case the kinematics of the low ion absorption profiles are 
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determined by the form of the rotation curve, V ro t(r), the 
value of the central column density, N±(0), and the thick- 
ness of the disks, h (see PW1). Because the lincwidth, Aw, 
increases with h we follow PW1 by adopting the largest 
plausible value, h — 0.3i?d, in order to maximize Aw. As 
stated in § 2.3, V rot (r) either equals V200 or /v^oo- We 
also make two assumptions about 7Vi(0). Either N±(0) 
equals the central column density in adiabatically con- 
tracted NFW halos, N$* MW , or N ± (0) = 10 212 ' cm" 2 
where the latter value is meant to illustrate scenarios in 
which star formation has consumed most of the disk gas 
(the value of Rd does not affect low ion kinematics since 
they are independent of the absolute scale length [PW1]). 
As a result the kinematics of the low ions are represented 
by four independent models. Because each of these will 
be linked to two kinematic models of the ionized gas (dis- 
cussed in § 5.2..1) we consider a total of eight kinematic 
models for each cosmogony. Their properties are summa- 
rized in Table 2. The models are designated by 4 letters; 
the first, M or N, specifies whether N ± (0) = N^ MW 
or 10 212 cm~ 2 , the second, V2 or fV, specifics whether 
V ro t(r) = V200 or /yV200! an d the third, R or B, indicates 
whether the velocity field at r < r cross is dominated by 
random motions or ballistic infall. Thus the model MV2B 
has a disk with MMW central column density, disk rota- 
tion speed given by V200! an d ballistic infall at r < r cross . 
Because the low ion kinematics of models MfVB, MV2B, 
NfVB, and NV2B are equivalent to those of models MfVR, 
MV2R, NfVR, and NV2R we shall discuss low ion results 
only for the former group. 

For a given V200 we find our single-disk CDM models 
with V ro t(r) = /VV200 to result in Aw that are larger than 
predicted by Kauffmann (1996) who did not correct for 
the gravitational contributions of the disk nor for adiabat- 
ically contracted halos to the disk rotation curves. Even 
so, none of our CDM models reproduces the observed low 
ion Av distribution. This is contrary to the expectations of 
MMW who conjectured that the higher V rot (r) produced 
by the more realistic rotation curves would result in a Av 
distribution compatible with observation. The reasons this 
does not occur are illustrated in Figure 2. In every case the 
median V200 of the the intercepted disks is much less than 
100 km s _1 . Because the Av predicted for rotating disks 
typically equals V200/3, the predicted median Av will be 
much less than 50 km s _1 . By contrast the median Av of 
the observed distribution is about 80 km s _1 . 

Figure 5 shows the results for model NfVB in each CDM 
cosmogony. With V ro t(r)=f v V 2 oo and N ± (0) = 10 212 cm~ 2 
this model yields the best case results because it generates 
the largest Aw of the 4 models. The larger velocity widths 
follow from (a) the higher rotation speeds, and (b) because 
the relatively low central column density of model NfVB 
results in smaller impact parameters out to the threshold 
column density iV(H I)=2xl0 20 cm -2 , and smaller impact 
parameters cause larger Aw (PW1). Nevertheless, appli- 
cation of the KS test shows that none of these models is 
likely to fit the data (see Figure 5). This conclusion also 
holds for ICDM even though ICDM predicts a larger frac- 
tion of high V200 halos than the other CDM models. This 
is because ICDM is also an hierarchical model (Peebles 
1999a, 1999b) and as a result a large fraction of the pro- 
togalactic mass distribution is in halos with V200 < 100 



km 



4.2. TF 



In the TF cosmogony we either let N± (0) equal the em- 
pirically determined Nq or or 10 212 cm~ 2 . As discussed 



V m 



in all TF models, where the 



previously V rot (r) 
distribution of V max is shown in Figure 2. Figure 5 also 
shows the model NfVB results for TF. Here the Aw test 
yields P^s(Aw) = 0.12. While lower than Pks(Aw) = 0.65 
exhibited by models in which every halo has V ro t(r) = 250 



km 



(PW1), this KS probability is sufficiently large 



that the more realistic TF model cannot be excluded. We 
shall check the robustness of these results for the TF and 
CDM models in § 6. 

5. RESULTS OF MONTECARLO SIMULATIONS: IONIZED 

GAS 

5.1. Normalization 

In order to fit the simulations to the data it is neces- 
sary to specify the function Nciv(b) where b is the im- 
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Fig. 5. — Comparison between empirical and Monte Carlo distri- 
butions of low ion test statistics for SCDM, ACDM, OCDM, ICDM, 
and TF cosmogonies. Results for NfVB model compare distributions 
of (a) Av, (b) {mm, (c) / e d 9 , and (d) ftyk t es t statistics which are 
defined in paper I. Briefly stated these quantities are the (a) ab- 
sorption velocity interval, (b) the normalized difference between the 
mean and median velocity, (c) the difference between the velocity of 
the strongest component and mean velocity, and (d) the difference 
between the velocity of the second strongest component and mean 
velocity. The quantity Pks is the KS probability that model and 
empirical distributions are drawn from the same parent population. 
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MODEL 

MfVB 

MV2B 

MfVR 

MV2R 

NfVB 

NV2B 

NfVR 

NV2R 



Table 2 
Disk-Halo Models 



Vrot = 

Ballistic 



N 



^200 

Random 



N MMW Qr N C or 



Vrot = /V V200 

Ballistic Random 



AT =10 21 - 2 cm- 


-2 


Vrot = V200 


Vrot=/v ^200 


Ballistic Random 


Ballistic Random 



"Central column density perpendicular to the disk given by eqs. (3) and (8) for CDM models (Nq IMW ) and by eq. (6) (N ~ !or ) for TF models 



pact parameter. This function is crucial as it constrains 
kinematic quantities such as Avcw f° r the C IV profiles 
by fixing the number of clouds per line of sight (see eq. 
30). We normalize Nciv{b) by determining the parame- 
ters K and (3 (eqs. 28 and 29) from comparisons between 
model and empirically determined frequency distributions 
of C IV column densities, f(Nciv, z )- The latter is the 
product of the number of damped Lya systems per unit 
"absorption distance", dn/dX, times g{Nciv), the con- 
ditional distribution of C IV column densities given the 
presence of a damped Lya system. We find that 



f(N, 



civ 



dX 

dz 



dn dX 
dz \ dz 

(1 + z) 2 



ANciv) 



[(l + z) 2 (l + n M z)-z(z + 2)Q A 



1/2 



(31) 



Because damped Lya systems are H I selected, the func- 
tion g will depend on the differential area of the inclined H 
I disks giving rise to damped Lya absorption. As a result 
the g will depend on impact parameter, 6, and hence on 
K and (3 through eqs. (28) and (29). We adjust K and (3 
by comparing model and empirical f(NciVi z)'s. 

We constructed the empirical g(Nciv) from the 32 C 
IV column densities inferred from the profiles in Figure 
1 in Paper I (The actual column densities are reported 
in Prochaska & Wolfe 1999 and in Lu ct al. 1996). To 
obtain the empirical f(Nciv, z ) we adopted the ACDM 
model and let dn/dz — 0.22 which is appropriate for the 
mean redshift of this sample. The results are shown as 
points with error bars in Figure 6. We compare this with 
the predictions for models MfVB and NfVB in the case of 
a ACDM cosmology. The results are valid for all MXXX 
and NXXX models respectively. This is because f(Nciv) 
depends on the distribution of impact parameters, but is 
independent of rotation speed and infall kinematics. In 
both classes of models, (3 = 1.5 provides a good fit to the 
data, while K = 2x 10 15 cm" 2 for model MfVBjmd 2x 10 14 
cm" 2 for model NfVB. When MM computed N C iv(b) for 
cool halo gas photoionized by background radiation, they 
found K and (3 to vary with V2007 i n contrast to our as- 
sumption of uniform K and (3 for halos of all mass. For- 
tunately, the C IV velocity profiles are independent of K 
because both Nqiv and N^jf v are linearly dependent on 
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Fig. 6. — Comparison between empirical and model CIV column 
density distributions. The points with error bars are data. The 
solid and short dashed curves are predictions for MfVB and NfVB 
models, /3 = 1.5 and the ACDM cosmogony. The short dashed curve 
corresponds to the MM model for halos with V200 > 150 km s — 1 
and in which (3 RJ —1 



K, and as a result K drops out of the determination of n c id 
which is crucial in determining the velocity profile widths 
(see eq. 30) . On the other hand the profiles do depend on 
(3. Figure 6 shows that (3 = — 0.9, which corresponds to 
the MM results for halos with V200 > 150 km s _1 , results 
in poor fits to the data. Therefore, the (3 s=s 1.5 assump- 
tion is consistent with the C IV data and is used in the 
calculations which follow. 

5.2. Test Statistics 

5.2.1. CDM Models 

In Figure 7a we compare the empirical Avqw distribu- 
tion with predictions by the ACDM cosmogony. We only 
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show results for the case V ro t(r)=V2oo, since halo kinemat- 
ics should be unaffected by disk rotation speed for halos 
of a given mass. We ignore results for the f m m, fedg, 
and ftpk statistics because in the case of C IV kinemat- 
ics Awciv is the most sensitive test statistic for testing 
any of the models. The principal difference between the 
distributions in Figure 7a stems from the different veloc- 
ity fields at r < r cross . Models with ballistic infall (MV2B 
and NV2B) predict larger median Awciv than models with 
random motions (MV2R and NV2R). For a given N±(0) 
the ballistic models are in slightly better agreement with 
the data because the larger Avciv's are closer to the ob- 
served values. 
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Fig. 7. — Comparison between empirical and Monte Carlo dis- 
tribution of A^civ- ( a ) Results for ACDM cosmogony and MV2B, 
MV2R, NV2B, and NV2R models, (b) Results for SCDM through 
ICDM cosmogonies and NfVR model, and for TF cosmogony with 
NV2R and NfVR models 



The value of iVj_(0) also causes differences between the 
distributions. Comparison between ballistic infall models 
MV2B and NV2B shows that MV2B predicts lower me- 
dian Auciv than NV2B. This is because the higher N±(0) 
of model MV2B results in larger impact parameters then 
cause the sightlines to sample the halos at larger radii 
where the infall velocities are smaller. In this case the 
lower At>civ of the MV2B model is in better agreement 
with the data. At the same time model MV2R is in better 
agreement with the data than NV2R, because the smaller 
impact parameters predicted by the latter model result 
in more sightlines traversing the r < r cross region where 
random motions give rise to Avcw that are also lower. 

The third difference between the Av distributions de- 
pends on the assumed cosmogony. In Figure 7b we use 
the NfVR model to illustrate the effects of the assumed 
cosmology. The cosmogonies differ according to the frac- 
tion of halos with large V200, a fraction which increases 



along the SCDM — > ICDM sequence. In every cosmogony 
in this Figure a "spike" in the Awciv distribution at A^civ 
= 50 km s _1 is present and increases in strength along 
the sequence. The spike arises from sightlines traversing 
the outer regions of high-mass halos at large impact pa- 
rameters. The high pressure of hot gas in massive halos 
(see Figure 3) compresses the clouds thereby reducing their 
cross-sections. For clouds of fixed mass the result is an in- 
crease in column density. In most cases only one cloud is 
required to satisfy eq. (30) at large b where Nciv is small. 
The spike occurs at Avcw = 50 km s _1 because this is 
the FWHM of the profile caused by a single cloud with 
an assumed internal velocity dispersion, o~i n t — 25 km s _1 
. In fact this value of <ii n t was chosen to reproduce the 
narrowest C IV profiles in our sample. 

It is worth noting that despite their differences most of 
the models yield large Pks(Avctv) values. This is in con- 
trast to tests of low ion kinematics. In that case all of the 
Pks(Av\ ow ) values were too small; i.e., none of the models 
was compatible with the observed distribution of low ion 
Aw. This could imply that while the infall interpretation 
for the high ions is correct, the disk interpretation for the 
low ions is incorrect. 

5.2.2. TF Models 

Figure 7 b also compares the data with the predictions 
of the TF cosmogony. In this case V200 equals V max or 
V-max/ fv where the V max are the input disk rotation speeds. 
As a result the CDM degeneracy of Avcw with respect to 
disk rotation speed is broken in the TF models. We show 
results for some examples to illustrate this effect. As ex- 
pected the Av — 50 km s _1 spike is highest for the NV2R 
model where V200 = V max ■ Because of the larger fraction 
of halos with high V200, the median Ati's are higher than in 
the CDM models. Agreement with the empirical Avcw 
distribution improves in the case of random motions at 
r < r cross and when V 2 oo = V max /fv] i-e., with the NfVR 
model (note the higher Pks(Avcw) value of the NV2R 
model is an artifact due to the large amplitude spike at 
Avcw = 50 km s _1 ), but the Pks(Avcw) values are 
still leass than 0.05. 

5.3. Correlation Tests 

In Figure 8 we compare empirical and predicted distri- 
butions of 5v and f ra uo- These are the differences between 
the mean velocities of the C IV and low ion velocity profiles 
and the ratio Avcw / Av\ ow , respectively . Comparison be- 
tween empirical and predicted cross-correlation functions 
for the C IV versus low ion velocity profiles is also shown. 

• Sv Test: 

In CDM all the models pass the 6v test at more than 
83% confidence. In Figure 8a we use ACDM to illustrate 
the effects of (a) disk rotation speed, (b) impact param- 
eter, and (c) halo velocity field. Comparison between 
the MfVB and MV2B models shows that disk rotation 
speed is the most important effect. Specifically Pks(3v) 
increases significantly when V200 decreases from /VV200 
to V2oo- This behavior is straightforward to explain. In 
the disk-halo models the half-width of the 6v distribution 
is roughly equal to the sample average of V ro t(r)sm(i). 
This is because the prototypical C IV profile comprises 
two widely separated absorption components symmetri- 
cally displaced about the systemic velocity of the halo. As 
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a result the mean velocity of the C IV profile equals the 
systemic velocity. On the other hand the low ion profile 
consists of multiple contiguous components comprising a 
single feature that is displaced to either side of the sys- 
temic velocity. In this case the mean velocity of the low 
ion profile is separated by ~ V ro t(r)sin(«) from the sys- 
temic velocity. Consequently the w 75 km s _1 half- width 
of the observed distribution limits the fraction of rapidly 
rotating disks. This constraint is especially severe for the 
TF models where rotation speeds exceeding 200 km s _1 
are typical. 
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Fig. 8. — Correlation Tests, (a) Comparison between empiri- 
cal and Monte Carlo distributions of 8v for ACDM cosmogony and 
MfVB, MfVR, MV2B,and NfVB models. The NfVB model is used 
with the TF cosmogony, (b) Same as (a) but f ra tio is substituted 
for 8v statistic, (c) Comparison between empirical and Monte Carlo 
predictions of t;(v) for C IV versus Low. Same models as in (a) and 
(b) are used. Curves are data (dark solid curve), ACDM MfVB 
model (dark dashed curve), ACDM MfVR model (dark dot-dashed 
curve, ACDM MV2B model (light solid curve), ACDM NfVB model 
(light dashed curve), and TF NfVB model (light dot-dashed curve) 

On the other hand the effects of impact parameter are 
not as significant. This is because impact parameter af- 
fects the width of the velocity features rather than the 
location of their velocity centroids. This explains why the 
the MfVB and NfVB results are so similar. Furthermore, 
the effect of halo velocity field is even less important as 
shown by comparison between the MfVB and MfVR re- 
sults. This tells us that for a sufficient number of C IV 
clouds, the location of the velocity ccntroid of the C IV 
profile is independent of whether the clouds are infalling 
or moving randomly. Therefore, Sv is set by the magnitude 
oiV rot (r). 

• J ratio J-CSt: 

A natural consequence of the disk-halo hypothesis is the 
prediction f ra tio > 1- Because the infall velocities of the 
high ions and the rotation speeds of the low ions occur in 



the same potential well, they both scale linearly with V2oo- 
However, owing to projection effects, the line-of-sight ve- 
locity gradients due to radial infall will exceed those due 
rotation. Thus, Auciv will be larger than A«i ow . But 
the ratios are too large in most CDM models because of 
the small A«i ow and the larger Avciv- As a result the 
best case models are those with large V ro t (r) for the disks 
and random motions at r < r cross for the halos. The best 
model is MfVR, as shown in Figure 8 b. The TF models 
also produce f ratio that are too high because in most halos 
the sightlines intercept the region r > r cross where large 
infall velocities are present. 

• Cross-Correlation Function: 

None of the models, neither CDM nor TF, predicts a 
cross-correlation function with large enough amplitude to 
fit the data. The reason is insufficient overlap in velocity 
space between the C IV and low ion profiles. In Figure 8 c 
the best case model is MfVR, indicating more overlap oc- 
curs when the C IV clouds experience random motions at 
f < r cross . Even better agreement is obtained with model 
MV2R (not shown) implying that overlap increases when 
the rotation speed of the disk is reduced. This interpreta- 
tion is supported by the results for the TF models which 
exhibit the worst agreement with the data. The TF mod- 
els predict the largest fraction of halos with high V ro t(r) 
and lowest fraction of low-mass halos in which random 
motions dominate the the velocity field at r < r cross . 

5.4. Model Summary 

In summary we find the following results for the models 
we haved tested so far: 

• While the extent of the low ion Av distribution rules 
out the single-disk semi-analytic CDM models, it is com- 
patible with the TF models (see also Jedamzik & Prochaska 
1998). 

• The C IV Av distribution is compatible with most of 
the CDM models. In the example shown in Figure 7 the 
best agreement with the data occurs for models with (a) 
ballistic infall at r < r cross and (b) central disk column 
densities given by N^ MW . The best agreement with the 
TF models is for ballistic infall, N =N^ MW , and V 2 oo = 

Vmax/fv- 

• The Sv distribution is compatible with all of the CDM 
models, but is too narrow for the TF models. This is a 
reflection of the low V ro t(r) predicted by CDM and the 
high V ro t(r) predicted by the TF models. 

• Most of the CDM models predict f ratio distributions 
with median values that are too large. This stems from 
the low values predicted for V ro t(r). The same problem 
holds for the TF models, but in this case the high f ra tio 
stems from the high values of Awciv- In both cosmogonies 
some models cannot be ruled out. 

• None of the models predicts C IV versus low ion cross- 
correlation functions in agreement with the data. 

6. PARAMETER TESTS 

The results of the last section may be summarized as 
follows: When the free parameters of a given model are 
adjusted to satisfy one test, the model inevitably fails a 
different test. If this is a generic feature of the disk-halo 
models, then they may not apply to the damped Lya sys- 
tems. This is an important conclusion and we wish to 
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determine how robust it actually is. Indeed the models 
are characterized by several free parameters that are not 
well determined, and it is possible we have not found the 
optimal set. For this reason we now investigate the sensi- 
tivity of our conclusions to variations of these parameters. 
In a series of trial runs we found the model kinematics 
to be most sensitive to three parameters. The first is the 
central perpendicular column density, iVj_(0). Because the 
range of impact parameters is limited by N±(0), it influ- 
ences both Ai>civ and Ai>i ow . To test the dependence of 
the kinematics on this quantity we simulate H I disks with 
iVj_(0) ranging between 10 208 and 10 222 cm~ 2 . These val- 
ues cover the low column densities of the NX XX models 
and approach the high values of the MX XX models. The 
second parameter is V cut , the low-end cutoff to the distri- 
bution of input V200 i n the CDM models or of input V max 
in the TF models. The kinematic results should be sen- 
sitive to V cu t especially in the case of CDM where V cu t is 
at the peak of the V200 distribution (see Figure 2). The 
value used in our models, V cu t = 30 km s _1 , is imposed 
by thermal expansion of photoionized gas out of the po- 
tential wells of halos with V200 lower than this (Thoul & 
Weinberg 1995). However feedback due to supernova ex- 
plosions might drive gas out of disks with V200 as large 
as 100 km s" 1 (see Dekel & Silk 1986, but see MacLow 
& Ferrara 1999). For these reasons we let V cu t vary be- 
tween 30 and 120 km s _1 . The third parameter is the C 
IV column density per cloud, N^f v . The definition given 
in eqs. (25) — (27) is for a spherical cloud of a given mass 
and C +3 /H ratio, i.e., Xqiv- To account for deviations 
from spherical symmetry or from our definition of Xciv 
(see eq. 26) we introduce the parameter qciv which is the 
ratio of the true C IV column density to our model defi- 
nition. We let qciv vary between 0.1 and 1.3. In order to 
supply the total C IV column density required at a given 
impact parameter, the number of clouds must increase as 
qciv decreases. This affects the C IV kinematics because 
Ai>crv will increase with cloud number. 

6.1. CDM 

The results of the parameter tests are summarized in 
Figure 9 which shows iso-probability contours in the V cu t 
versus N±(0) plane. The contours correspond to 0.01, 
0.05, and 0.32 for P KS (Av low ) (9a), P KS (Av CIV ) (9b), 
Pks(Sv) (9c), and Pusifratio) (9d). We show results for 
qciv = 1-0 since smaller values are found to result in f ra tio 
that are too large. We choose a variant of the NfVR model 
in which -/Vj_(0) is a free parameter. We abandon ballistic 
infall in favor of random motions at r < r cross because we 
find that random motions produce more overlap in veloc- 
ity space between low ion and C IV profiles and as a result 
they produce better agreement between model predictions 
and the empirical test statistics. 

Figure 9a shows the results for the low ion Av test. As 
expected models with the standard value V cu t = 30 km s _1 
are improbable for reasonable values of N±(0). Rather 
V cu t increases along the Pjfs(Avi ow ) = 0.05 contour from 
35 km s" 1 at logiVx(O) = 20.8 cm~ 2 to 105 km s" 1 
at logiVx(O) = 22.0 cm~ 2 . An increase in N±(0) means 
larger impact parameters which in turn imply smaller Av\ ow 
(see § 5.1). Therefore, an increase in V cu t must accompany 
the increase in N±(0) to boost the fraction of high-T^oo ha- 
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Fig. 9. — Iso-probability contours in V cu t versus N±(Q) plane 
resulting from parameter search of ACDM cosmogony where NfVR 
model is assumed in which is a free parameter and where qciv = 
1.0. Contours correspond to Pks = 0-32 (solid lines), 0.05 (dots), 
and 0.01 (dashes). Tests are (a) Ai>i ow , (b) A^ciVi ( c ) Sv, and (d) 

Jratio- 



los required to maintain the extent of the A-ui ow distribu- 
tion. By contrast Figure 9b shows that V cu t hardly varies 
with N±(0) along the Pks(Avciv) contours. In our CDM 
models most sightlincs traverse halos with low V200 where 
cloud motions in NfVR models are dominated by random 
velocities drawn from a Gaussian with dispersion that is 
an insensitive function of radius (eq. D4). Consequently 
Awciv will be independent of impact parameter and hence 
independent of iVj_(0). Therefore, V cu t need not vary with 
N± (0) to maintain the extent of the Avqw distribution. 

Figure 9c and d shows the iso-probability contours for 
the 5v and fratio tests. In Figure 9c none of the con- 
tours rises to Pks{5v) — 0.32 in the V cu t versus N±(0) 
plane. The shape of the Pks{5v) = 0.05 contour has the 
following implications. First, the insensitivity of the con- 
tour to N±(0) just indicates that the displacement of the 
low ion velocity centroid from the systemic velocity of the 
galaxy is determined by rotation speed rather than At>i ow . 
Second, models with V cu t > 55 km s _1 are highly un- 
likely as they produce Sv that are too large. In Figure 
9d, -/Vl(O) decreases with increasing V cu t along all iso- 
probability contours. Because Auciv increases with in- 
creasing Vcut, Af_i_(0) must show a corresponding decrease 
to boost A«i ow . Otherwise f ra tio becomes larger than ob- 
served. 

Figure 10 shows the Pks — 0.05 contours from the 
previous Figure. The horizontal lines trace out the region 
in the V cu t versus N±(0) plane in which Pks > 0.05 for 
all 4 tests; i.e., the parameter space of acceptable mod- 
els. Physically, the resulting range of V cu t (i-e- 35 to 
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50 km s -1 ) is acceptable for models in which gas pho- 
toionized by ionizing background radiation escapes from 
low- mass halos (Thoul & Weinberg 1995; Navarro & Stein- 
metz 1997). On the other hand the upper limit on N±(0) 
(i.e., \ogN±(0) < 21.2 cm~ 2 ) may be too low to explain 
the shape of the H I column-density distribution function 
(see § 7). Furthermore, this Figure indicates these mod- 
els may not be viable as they occupy a small fraction of 
the depicted parameter space. Notice that the horizontal 
Sv contour at V cu t ~ 50 km s _1 is crucial in restricting 
the acceptable region to such a small area. Because Sv 
is set by rotation speed we investigated the NV2R mod- 
els to determine whether the lower rotation speeds would 
enlarge the acceptable region. The results are shown in 
Figure 1 1 . As predicted the lower rotation speeds lift the 
restricting Sv contour from ss 50 to w 95 km s _1 . How- 
ever, the lower rotation speeds also increase f ratio, with 
a consequent lowering of the Pr s(f ratio) contour in Fig- 
ure lid. Consequently there is no region in the V cu t verus 
N±(0) plane in which all 4 tests result in Pks > 0.05 for 
the NV2R model. Therefore, our conclusion concerning 
the size of acceptable regions in parameter space appears 
to be robust. 
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Fig. 10. — solid lines are contours corresponding to Pks = 0.05 in 
Figure 9. The region in which Pks > 0.095 for all 4 tests is denoted 
by horizontal lines 

Turning to the cross-correlation function we find that 
none of the models within the range of V cu t and -/Vj_(0) 
depicted in Figure 9 results in C IV versus low ion cross- 
correlation functions with acceptable \ 2 values. Appar- 
ently the combination of radial infall and disk rotation 
produce C IV and low ion absorption profiles with insuffi- 
cient overlap in velocity space to explain the data. Better 
agreement is obtained when logiV^(O) > 23 cm -2 . But 
column densities this high are ruled out by the other tests. 
We shall return to this dilemma in § 7. 



Fig. 11. — Same as Figure 9, except with NV2R model 



In summary, by varying V cu t, JVj_(0), and qciv we find 
regions of parameter space where the ACDM models are 
in better agreement with the data than for the "standard" 
values of the parameters adopted above. This is especially 
true for the low ion Av test which ruled out most of the 
"standard" models (where V cu t = 30 km s _1 ). However, 
the models may still not be viable because of the restricted 
range of allowable parameters, and because none of the 
models is compatible with the C IV versus low ion cross- 
correlation function. 

6.2. TF 

The corresponding results for the TF models are shown 
in Figure 12. In this case the best fit value of qciv equals 
1.3. Figure 12a illustrates the results for the low ion Av 
test. The Figure shows the Pks{^Iow) — 0.05 contour 
to enclose a larger area of parameter space than in the 
CDM case. As in CDM, V cu t increases with N±(0) along 
iso-probability contours. However, the A«i ow contours are 
less sensitive to V cu t because the input halo distribution 
does not peak at V cut as in CDM (see Figure 2). Figure 
12b shows that by contrast with CDM the Pks(Avciv) 
contours are sensitive functions of N±_ (0) . This is because 
in the TF model more sightlines traverse high V200 halos 
where C IV clouds undergo infall at r > r cross , and as a 
result Awciv lS a sensitive function of impact parameter, 
and therefore of N±(0). 

The results for the Sv and f ratio tests are plotted in Fig- 
ure 12c and d. All of the models result in Pks{5v) < 0.05. 
Clearly the large fraction of rapidly rotating disks encoun- 
tered in the TF model produces Sv that are too high. By 
contrast the results for the f ratio lead to Pks(Sv) > 0.05 
throughout the parameter space depicted in the Figure 
except for the upper right portion where Pks(8v) > 0.32. 
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Fig. 12. — Same as Figure 11, but for TF cosmogony 



But part of the improvement here is caused by spikes in the 
Awciv distribution near 50 km s" 1 which prevent f ra tio 
from exceeding observed values when the At>i ow are small. 
While the physical basis for such spikes is understood (see 
§ 5.2.1), they have not been confirmed by the data. As in 
CDM, none of the TF models produce C IV versus low ion 
cross-correlation functions with acceptable \ 2 values. In 
fact the higher fraction of disks with large rotation speeds 
predicted by the TF model increases the displacement be- 
tween low ion and C IV profiles which produces even lower 
cross-correlation amplitudes than in the CDM models. 

We also tested the sensitivity of the TF model results 
to variations of (a) 7 and V* , the power-law exponent and 
fiducial rotation speed in the Tully-Fisher relation (see eq. 
13), and (b) the Schecter function exponent, a (see eq. 
15). Within the parameter range, 3 < 7 < 4, (cf. Gio- 
vanelli et al. 1997; Sakai et al. 2000) the test statistic 
distributions do not change significantly. By contrast, in- 
creasing a above 1.2 reduces At>i ow significantly owing to 
the larger fraction of low mass halos. This would occur if 
the damped Lya galaxies were drawn from the luminosity 
function measured for the Lyman-break galaxies (Stcidel 
et al. 1999), since in that case a = 1.6 (see § 2.1). But 
overall the improvement of the model is poor because f ratio 
increases to unacceptably high values. The results also 
change when we vary V* above 280 km s _1 or below 220 
km s _1 . But in both cases the model does better against 
some tests and worse against others. In no case did the TF 
model pass all 4 tests at more than 95 % nor did the fits of 
the cross-correlation function become acceptable. In that 
sense the TF models are in worse agreement with the data 
than the CDM models. 

7. SUMMARY AND CONCLUDING REMARKS 



We used accurate kinematic data acquired for a sample 
of 35 damped Lya systems to test the standard paradigm 
of galaxy formation; i.e., the scenario in which galaxies 
evolve from the dissipative collapse of virialized gaseous 
halos onto rotating disks. The data was presented in the 
form of velocity profiles of high ions, intermediate ions, 
and low ions in Paper I. In this paper we considered semi- 
analytic models, specifically the MMW models in which 
ccntrifugally supported exponential disks are located at 
the centers of dark matter halos drawn from mass dis- 
tributions predicted by standard CDM cosmogonies, and 
where infall of ionized gas from the halo occurs. We also 
considered the null hypothesis that current disk galaxies 
were in place at z > 3 (the TF models). We tested the 
models with Monte Carlo techniques by comparing dis- 
tributions of test statistics generated from observed and 
model velocity profiles. We utilized eight test statistics: 
A^iow, fedg, fmm, and f tpk for the low ion profiles (see 
Paper I), Awqiv for the C IV profiles, and 5v, f ratio, and 
£(v) for comparing low ion and C IV profiles. 

First, we discuss the general implications of our work. 
As discussed in Paper I velocity profiles overlapping in 
velocity space in such a way that Avcw > A«i ow are nat- 
urally reproduced by scenarios where low ion and high ion 
kinematic subsystems are in the same gravitational po- 
tential well. In the collapse scenario the velocity fields of 
both subsystems scale as V200, yet more of V200 is pro- 
jected along the line of sight by gas undergoing radial 
infall than by gas confined to rotating disks. This was 
confirmed by our Monte Carlo simulations of radial infall 
of ionized gas clouds onto neutral rotating disks. Indeed 
in some cases the infall velocities exceed V200 resulting in 
A^civ/Awiow ratios that are too large. By contrast, sce- 
narios in which the high ions are embedded in gaseous out- 
flows (e.g. Nulsen et al. 1998) or any flows not generated 
by dark-matter potentials determining low ion velocities 
will in general not satisfy these constraints. 

Next we discuss specific conclusions arising from this 
work, in particular the results of model testing. Tests of 
models with the standard parameters discussed in § 4 and 
§ 5 led to the following conclusions: 

(1) In the case of the low ion gas none of measured dis- 
tributions of Aw low , f edg , f mm , and f tpk were compatible 
with the predictions of the CDM cosmogonies at the 95% 
confidence level. By comparison, the TF model was com- 
patible with the data at 88 % confidence level. 

(2) For the high ion gas we considered only the Awciv 
distribution. Comparison with the data showed CDM 
models with the high column densities predicted by MMW, 
i.e., iVj_(0) = N±(0) MMW , were in good agreement with 
the data. CDM models with significantly lower N±(Q) 
were not as good because they produced overly large A^civ- 
For the same reasons TF models with high -/Vj_(0) were in 
better agreement with the data than with low N±(0). 

(3) To test model predictions for the relative properties 
of the high ion and low ion gas we considered the 5v distri- 
bution. The CDM models were compatible with the data 
while TF models were not. Apparently the high rotation- 
speed TF disks displace the asymmetric low ion profiles 
too far from the velocity centroids of the C IV profiles. 

(4) Tests of the f ratio distributions showed neither CDM 
nor TF model predictions agreed with the data at 95 % 
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confidence. 

(5) Neither the CDM nor TF models predicted C IV 
versus low ion cross-correlation functions that were com- 
patible with the data at 99 % confidence. 

We then explored parameter space to determine whether 
these conclusions were robust (see § 6). We varied three 
crucial parameters: V cu t, ATj_(0), and qciv- This exercise 
led to the following conclusions: 

(6) Figure 9 shows that the NfVR-ACDM model is com- 
patible with the A«i ow , Auciv, Sv, and f ra uo tests at more 
than 95 % confidence throughout the small area of the V cu t 
versus N±(0) plane shown in Figure 10. This is a serious 
shortcoming since the NfVR model used for the compar- 
ison is "optimistic" in that it predicts a constant veloc- 
ity rotation curve with maximum rotation speed V ro t(r) 
= ,/VV2oo- It also predicts low impact parameters owing 
to the restriction iVi(O) < 10 212 cm~ 2 . In addition, a 
small decrease in disk thickness would eliminate consis- 
tency throughout parameter space; i.e., the disks must be 
thick. 

(7) Figure 12 shows the NfVR-TF model to be compati- 
ble with the At>i ow , Avcw, and f ra uo tests and incompat- 
ible with the Sv test at the 95 % confidence. 

(8) Neither CDM nor TF models produce C IV versus 
low ion cross correlation functions that were consistent 
with the data in the parameter space shown in these Fig- 
ures. 

What have we learned from the model tests? Because 
the CDM models pass 4 out of 5 tests and the TF models 
pass 3 out of 5 tests, the CDM models appear to be more 
plausible. But to achieve this result it was necessary to 
adopt a flat rotation curve with V ro t{r) = /yV^oo- This 
is the maximum rotation speed possible for a model disk, 
and V ro t(r) in realistic protogalactic disks are probably 
lower. However, Figure 11 shows that a parameter search 
for models with V ro t(r)=V2oo reveals no region in parame- 
ter space which is compatible with the 4 kinematic tests at 
95 % confidence. Second, the limit N±(0) < 10 212 cm~ 2 
indicates that exponential disk models should predict a 
steepening of the column-density distribution function at 
N > 10 212 cm~ 2 . This effect is not present in the data 
(Wolfe et al. 1995; Rao & Turnshek 1999). Third, the 
failure of any model to reproduce the C IV versus low ion 
£(v) indicates significant overlap in velocity space between 
the low ion and high ion velocity fields was not achieved. 
Fourth, the CDM models predict that most of the damped 
Lya systems occur in low mass halos where the kinematic 
state of the ionized gas is highly uncertain (see § 3.1). 

Does this mean that disk-halo models for damped Lya 
systems are ruled out? We think it is premature to reach 
this conclusion. Rather we take these results to mean that 
if the collapse scenario is correct, a stronger coupling be- 
tween the kinematic subsystems is required. One possi- 
bility that comes to mind is for low ions to be associated 
with the infalling C IV clouds. This would increase the 
low ion line-of-sight velocities and cause smaller differences 
between the C IV and low ion velocity profiles. But the 
problem is there is no evidence for low ions with high ion 
kinematics. Another way to couple the subsystems is to 
include the angular momentum of the halo gas. This is ne- 
glected in the radial infall model. As the C IV clouds ap- 
proach the disk they spin up and experience azimuthal ve- 



locity components approaching V ro t(r). The idea is plau- 
sible if the angular momentum vector of the infalling gas 
is related to that of the disk, and if clouds near the disk 
are likely to be detected. It is encouraging that recent N- 
body simulations show the angular momenta of disk and 
halo to be correlated (Weinberg 2000). It is also encour- 
aging that the density of clouds along the line of sight is 
highest near the disk. Still, if this idea doesn't work one 
would be forced to abandon the disk-halo hypothesis; i.e., 
one of the standard paradigms of galaxy formation. 

Can the kinematic data be better explained by scenarios 
other than infall of ionized gas onto rotating disks of neu- 
tral gas? First, we already discussed problems associated 
with outflow models. Second, lacking analytic expressions 
for the various quantities, it is not clear whether numerical 
simulations of damped Lya systems are compatible with 
the kinematics of the ionized gas. However, the density 
contours in Haehnelt et al. (1998) show the C IV clouds 
to be within ~ 10 kpc of the low ion clouds in which case 
both would be subjected to the same dark-matter gravi- 
tational field. Because the low ion gas is not confined to 
rotating disks, it is not obvious why the predicted J 'ratio 
should exceed the observed lower limit of unity. To satisfy 
this constraint one must consider contributions to the C IV 
profiles from gas outside the dark-matter halos, perhaps in 
the fashion described by Rauch et al. (1997) for the C IV 
QSO absorption lines. Third, the scenario described by 
the semi-analytic modeling of Mailer et al. (2000) may 
provide a good fit to the C IV kinematics. In particu- 
lar, both the low ion and C IV kinematics arise from the 
orbital motions of mini-halos accreted onto more massive 
halos with V200 ~ 150 km s _1 , implying Awciv > A^i ow 
and that C IV and low ion profiles are well correlated. The 
current difficulty with the model is to physically motivate 
the very large Mestel disks required to explain the low ion 
kinematics. In any case performing the tests outlined in 
this paper will reveal how robust these models actually 
are. 

We wish to thank the following people for valuable dis- 
cussions: Robert Benjamin, James Bullock, Eric Gawiser, 
Chris McKee, Jim Peebles, and Amos Yahil. JXP acknowl- 
edges support from a Carnegie Postdoctoral Fellowship. 
AMW was partially supported by NSF grant AST0071257. 

APPENDIX 

PRESS-SCHECTER THEORY FOR ICDM 
COSMOGONIES 

According to Press-Schecter theory the density of bound 
halos is in the mass interval (M, M + dM) is given by 



n(M, z)dM 



-2^dM 



(Al) 



MdM 

where p is the mean density of matter and F is the fraction 
of objects with masses < M that collapsed by redshift, z. 
For Gaussian distributed density fields, F is given by 



F = (1/2) 



1 - erf 



(A2) 



,x/2£>(X)Ao(M), 

where D(z) is the density-contrast growth factor (Peebles 
1980) and Aq(M) is the current rms density contrast with 
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mass M. While eq. (A2) is applicable to ACDM mod- 
els, it does not apply to ICDM models where the density 
contrasts are not Gaussian distributed. Rather, we shall 
use an analytic fit to results of numerical computations 
(Peebles 1998b) which shows that in ICDM, F is given by 



F = 0.37exp 



0.67 



5 C 



(A3) 
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For power spectra in which P{k) oc k n , Ao(M) ex M~(~e~ 
As a result 
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where M#, the mass going non- linear at redshift z, is ob- 
tained from the expression, A (M,) = S c /D(z) (see White 
1997). One finds that 



M,(z) = M s 



A (M 8 



D(z) 



(A5) 



where the mass in a sphere with radius R = 8/i _1 Mpc, 
Ms = 5.9 x 10 14 QmMq. Throughout the paper, we de- 
note Ao(Mg) by the more conventional symbol, as- In the 
case of Gaussian linear density fields applicable to ACDM 
models wc have 



ACDM : n{M, z)dM 



2 p 



dA a 



x exp 



7T M DA 2 dM 
S 2 



2D 2 A 2 



dM 



(A6) 



where in this case Ao(M) is determined by the full CDM 
power spectrum (Bardeen et al. 1986). 

HOT GAS IN NFW HALOS 

Following MM we assume the hot gas is in hydrostatic 
equilibrium with the dark matter potential and exhibits 
adiabatic temperature and density profiles out to r m i n = 
min(r coo ;,r2oo)- For the NFW rotation curves given in eq. 
(7) the density and temperature profiles of hot gas within 
r min are given by 



Ph(r) = Ph(r m in)[l + 0.8Z NFW (r)} 
T h (r) = T h (r mln )[l + 0.8Z NFW (r)} 



3/2 



where 



ZNFw(r) = 



Jnfw 



ln(l + x) hl(l + X m in) 



and 



Snfw = [ln(l + c) - c/(l + c) 



(Bl) 

; x = r/r s 
(B2) 

(B3) 



At r > r m i n the hot gas either follows isothermal profiles 
when r2oo > r coo i or does not exist when r2oo < f coo i . The 
temperature and density of hot gas at r m i n are given by 



-*- hymin ) -*-v 



and 



Ph{1 



fgV, 



g K 200 



AirGr 2 5 NFW x coo i{l + x cool ) 2 



(B4) 



i *Ecool — Tcoolf'Ts 



where f g is the mass fraction in gas. To solve for r coo \ we 
follow MM and let p h {r min ) = 5(pm H ) 2 V 2 oa 2 /4A(T mr )t M 
where A(T) is the cooling function. From the last equation 
we find that x coo i is the root of the following polynomial: 

Xcool(l + Xcooi) 2 - {c/5 NFW )(Rcooi/r s ) 2 = , (B6) 

where R CO oU the cooling radius of a singular isothermal 
sphere, is given by 



i^cool — 



A(T v )t M m d 



1/2 



(B7) 



hTiGp 2 m 2 H J 
(see MM). The root of interest is given by x CO oi > 1- 

COOL GAS 

The cool gas in the halo comprises photoionized clouds 
formed from hot phase gas at r < r m i n where r TO j„ = 
min(r coo ;,r2oo)- Therefore, the mass of gas in the cool 
phase is given by 



Mcool = .f g M{r min ) - M h (r min ) , 



(CI) 



where M(r m i n ), the total mass within r m i n , and M^(r TO j„), 
the mass of hot gas within r m i n , are given by 



and 



M(r n 



M h (r mm ) 



Train * rot vmin) 

G 



4nr ph{r)dr 



(C2) 



(C3) 



For NFW halos we use the rotation curve in eq. (7) 
to relate V ro t(r m in) to V200, cqs. Bl— B4 to compute the 
integral in eq. C3, and finally eq. B5 to evaluate Ph{fmin)- 
The result is 



Mcool = 



f g r 200 V200 



m^l T" Xjnirfij 



GSnfw 

iNFw(Xmin) 



[Xcool )\± ~r Xcool) 



(J- 1 Xmin ) 

;x = r/r s , (C4) 



where 

iNFw(Xmin) 



x 2 [l + 0.8Z NFW (x)]dx . (C5) 



MM assume the following form for the density of the cool 
gas: 
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Pc(r) = 



M, 



cool 



(C6) 



47rr 2 KtM ' 

where V c is a characteristic infall velocity. Combining the 
last equation with p c (r) — Ci/r 2 , we have 



d = 



,/>200^ 



2 
200 



AnGVctMSNFW 



ln(l + a; mi „) - 



^J- T" %>min) 



I NFW (x r , 



\Xcool)\^ T %cool) 



.(C7) 



Having obtained Ci we compute r c 
these results with eq. (17). We find r cr 
of the following equation: 



s by combining 
to be the roots 



r 2 [l + ZAT FW (r)] 



5/2 _ 



CaT c 







(C8) 



/^/i ymin )-L h \Tmin ) 

We compute the cloud cross-section, A(r), by combining 
these results with eq. 18. We hnd that 



In the case where the kinematics are dominated by random 
motions we solve the Jeans equations for a r (r), the radial 
velocity dispersion of the clouds, assuming locally isotropic 
velocity dispersions (Binncy & Tremaine 1987). In this 
case we have 



Orijf 



1 



v{r) J r 



dxv(x)g(x) 



(D3) 



where v(r) is the average density of clouds. Assuming v(r) 
ex p c (r) we find that 



o*(r) = V, 



2 I 

200 I 



;-!<) 



NFW 



2z 2 ln(a 



1 



- 2x 2 ln(l + x) + 



2x-\ 



= r/r 



,(D4) 



A{r) = 1T [ C 2 M 2 /|(r)[l + Z NFW {r)]-*/ 3 ;r > r. 
I C 3 M c 2 / d 3 (r cross )(r/r mi „) 4 / 3 ]r<r cro 



cross 
ss J 



where 



(C9) 



C 2 = 



3T C 



-i 2/3 



,c 3 = 



Q r 2 



-i 2/3 



47rCi 



47rp/i (r mm ) T/j (r min ) 

(CIO) 
We compute the cloud mass by integrating eq. (19) from 
Train to r to find M c i d (r) at r > r cross , and assume 
M c id(r) = M cld (r cross ) at r < r cross . We find that 



M cld (r) 



M c i d (r min ) 
M cid {r cross) ;r <r, 



1 + Ittjj^ i^: 7 " t 1 + 0.8^ Flv ( a ;)]- 1 /6 dx 



1 Jr J r 



:r > r c 



(Cll) 



CLOUD KINEMATICS 



The infall velocities of the clouds are computed as fol- 
lows. At r > r cross we compute V r (r) from eq. (22) by 
using the solutions for M c i d {r) given in Appendix C and 
noting that the gravitational acceleration is given by g(r) 
= V 2 ot (r)/r. From eq. (7) we find 



l(r) = K, 2 oo 



1 



1 



r s Snfw x 



ln(l + x)- 



l + x 



x = r/r s 



(Dl) 

At r < r cr0 ss we consider two scenarios. In the case 
of ballistic infall we assume M c i d (r) = const. For clouds 
with infall velocity V r (r cross ) at r = r cr0 ss, the ballistic 
solution for NFW halos is given by 



V™{x) 



V r V cross) + V200 



111^1 -\~ Across) 



2c 



Snfw 



ln(l + x) 



x = r/r s .(D2) 
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